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We explore the main processes involved in the evolution of gen- 
eral quantum systems by means of Diagrams of States, a novel 
method to graphically represent and analyze how quantum infor- 
mation is elaborated during computations performed by quantum 
circuits. We present quantum diagrams of states for representa- 
tions of quantum states by density matrices, partial trace oper- 
ations, density matrix purification and time-evolution by Kraus 
operators. Following these representations, we describe by dia- 
grams of states the most general transformations related to single- 
qubit decoherence and errors. 

Diagrams of states prove to be a useful approach to analyze quan- 
tum computations, by offering an intuitive graphic representa- 
tion of the processing of quantum information. They also help 
in conceiving novel quantum computations, from describing the 
desired information processing to deriving the final implementa- 
tion by quantum gate arrays. 
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1 EVOLUTION OF QUANTUM SYSTEMS BY GRAPHIC REPRE- 
SENTATION OF STATES 

We explore the main processes involved in the evolution of general quantum 
systems by means of Diagrams of States. The representation by diagrams of 
states is a novel method to graphically represent and analyze how quantum in- 
formation is elaborated during computations performed by quantum circuits. 
In the widely-used representation by quantum circuits, horizontal lines rep- 
resent single qubits constituting the considered quantum system. In contrast, 
in diagrams of states we draw a horizontal line for each state of the compu- 
tational basis. Therefore, diagrams of states are less synthetic in respect to 
quantum circuits, but allow a clear and straightforward visualization of the 
quantum information processing. 

We previously introduced this method by defining basic representations 
for standard quantum operations and providing examples of basic practical 
quantum computations |T|. In this paper, diagrams of states are applied to ex- 
plore the density matrix representation of quantum states, the time-evolution 
of open quantum systems by Kraus operators, and a general model for single- 
qubit decoherence and errors. As in previous related works lUEl, diagrams 
of states will be used both as a novel approach to investigate quantum compu- 
tations, in addition to (or in substitution of) standard methods like analytical 
study and Feynman diagrams and as an auxiliary tool to construct novel 
quantum computations from the desired manipulation of quantum states. 

This paper is organized as follows. In Section |2l we consider the repre- 
sentation of quantum states by density matrices, and their partial trace and 
purification processes for composite systems; as practical examples, we illus- 
trate physically feasible procedures to reduce density matrices in two-qubit 
systems and to purify density matrices of single-qubit systems. We conclude 
this section with the Kraus representation for the most general time-evolution 
of density matrices. In Section[3l we explore the most general transformation 
to model single-qubit decoherence and eiTors fSl. Finally, in Section |4] we 
present our conclusions. 

Throughout this paper, in order to perform the analysis of given quantum 
processes, we shall directly derive diagrams of states from the quantum cir- 
cuits associated with the physical implementation of the processes. These 
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diagrams can easily be rearranged into new simpler diagrams, which better 
visualize the overall manipulation of information from input to output: We 
shall refer to the former as complete diagrams and to the latter as simplified 
diagrams. 

Any sequence of logic gates must be read from left (input) to right (output), 
both for conventional quantum circuits and for their representations by means 
of diagrams of states. From top to bottom, qubits run from the least significant 
(lsb) to the most significant (MSB). 

2 DENSITY MATRICES AND THE KRAUS REPRESENTATION 

The density matrix representation (see, e.g., ||6|, pages 258-264, or Q, pages 
295-307) is very useful whenever representing non-pure states or describing 
their time evolution. Any quantum state can be represented by means of a 
density matrix which, in turn, can be expressed by its spectral decomposition: 



where are the density matrix eigenvectors and {A;} are the correspond- 
ing eigenvalues. 

2.1 Partial trace and reduced density matrices 

The Hilbert space of a composite system which is constituted by two inde- 
pendent or interacting subsystems can be obtained as the tensor product of 
the Hilbert spaces corresponding to the two constitutive subsystems. 

Given the density matrix of a general state of the composite system, a par- 
tial trace operation consists in finding a description of states which considers 
only elements related to a certain subsystem (denoted as subsystem /) while 
disregarding any element related to the remaining subsystem (denoted by //): 



where i is the index corresponding to subsystem II. 
Reduced density matrices in the spaces of two qubits 

For the sake of simplicity, we consider a system composed of two single- 
qubit subsystems, whose general overall state is described by a density matrix 
of dimension 4x4. By performing on this density matrix the partial trace 
operations on each one of the two single-qubit subsystems, it is possible to 
define the two corresponding reduced density matrices (of dimension 2 x 2). 




(1) 




(2) 
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The two possible partial trace operations are graphically represented in 
Figure [T] The processing of information involved in the partial trace oper- 
ations is most clearly illustrated in the diagrams of states (left). Evidently, 
the diagram related to the partial trace on the least significant qubit (lower) 
can be directly derived from the diagram related to the partial trace on the 
most significant qubit (upper) by appropriately permuting the qubits consti- 
tuting the system; the permutation is obtained by simply applying swap gates 
before and after the partial trace operation. By appropriate immersions and 
permutations of diagrams of states |1|, the method can be easily general- 
ized to describe partial trace operations in systems of higher dimensionality 
and in respect to subsystems constituted by any combination of qubits. Fi- 
nally, the diagram-of-states representation for partial trace operations will be 
a fundamental tool for the description of decoherence and evolution of open 
single-qubit systems, addressed in Section[3j 
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Figure 1 

Graphic representations of partial trace operations for a two-qubit system, on the most 
significant qubit (upper figures) and on the least significant qubit (lower figures). From 
left to right: quantum circuits, additions of matrix entries, diagrams of states. The 
reduced density matrix of the subsystem constituted by the least significant qubit is 
obtained by adding the 2x2 sub-matrices highlighted in gray shading. Similarly, the 
reduced density matrix of the subsystem constituted by the most significant qubit is 
obtained by adding in pairs the entries highlighted in gray shading. 



2.2 Purification of a density matrix 

The purification of density matrices (see, e.g., [6], pages 267-270) can be 
considered as a sort of inverse process to partial trace operations. A general 
non-pure quantum state in a system denoted by / is represented by its density 
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matrix . This density matrix can be seen as a reduced density matrix 
obtained by partial trace from a larger density matrix p, describing a pure 
state in a larger system. The larger system would include the initial system 
/ and a second ancillary quantum system //, and it is denoted by / + //. 
The purification of the density matrix p^ consists in determining the quantum 
system I + II which leads to the reduction process. 

By using the Latin set of indices ...} when referring to the initial 
quantum system / and the Greek set of indices {a, (3, ...} when referring to 
the ancillary quantum system //, a general pure state in the overall system 
I + II is described by: 

1$) =^ a„|z)|a) (3) 

i cx 

and by the density matrix: 

ia j 13 

The assigned density matrix p^ is in general not diagonal: 

k I 

By imposing that the density matrix p^ is obtained by partial trace of the 
density matrix p of the overall system in respect to the basis {I7)} of the 
subsystem //, we obtain: 

kl 7 [ i j J 

that is, a set of equations for the entries of the density matrix p^: 

Pij — ^ CijC*^, (7) 

7 

where the coefficients pjj are known, while the coefficients Ci-y are to be 
determined. 

As is well known, equations (Q are solvable if the space of the basis {I7)} 
has a sufficiently high dimensionality. That is, given any density matrix in a 
defined quantum system, it is always possible to consider an extended quan- 
tum system in which the assigned non-pure state can be obtained by partial 
trace of a pure state in the extended system. 
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Density matrix purification for a single-qubit system 

For the sake of simplicity, we consider the purification process for a density 
matrix p assigned in a single-qubit system. A single ancillary qubit is suffi- 
cient to obtain the corresponding pure state \^) in a larger quantum system. 
From equations (|7]) we obtain a three-equation system in four unknown 

Cia - 

Pll — CooCqo + CqiCqi 

Pi2 = CooC*Q + CoiCii ~ (8) 
P22 — CiqC^q + C'liC'li 
We can assume the parameter Cqo to be real, since a common phase is arbi- 
trary, and we can impose the condition Coi ~ 0, thus obtaining: 



n - /t; n (\ n P12 n / pn P22 - |pi2p 

\/Pll V 

(9) 

Thus, the following pure two-qubit state performs a purification of the given 
density matrix p: 



I*) = yPii I0)|0) + ^ 11)10) + /^li^Hi^kli^ (10) 

VPii V 

Figure |2] shows the quantum circuit to generate the quantum state for the 
purification of the density matrix p. The processing of information is clearly 
illustrated by the complete and simplified diagrams of states in figure [3] In 
these diagrams, the active information contained in the initial states is ma- 
nipulated from left to right along thick lines, while thin lines correspond to 
absence of information. 

After the application of the "A" gate array of Figures |2] and [3] we obtain 
the state: 

cos 6*1 


cos(?2 sin 6*1 
sin^?2 sin 01 

Purification is then achieved by applying the "B" gate array, where the con- 
trolled -phase gate has phase — arg(pi2). Thus, the pure state \^) of 
equation (fTOl i is generated by two-qubit unitary operations. 

2.3 The Kraus representation 

The most general time-evolution of a density matrix can be described by the 
Kraus representation ([8 1 or see, e.g., JS], pages 278-284). This representation 
is schematically illustrated in Figure |4] 



I*) = 



(11) 
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Figure 2 

Quantum circuit representing the syntliesis of a two-qubit state purifying the density 
matrix of a single-qubit system, by unitary operations performed on a two-qubit reg- 
ister. The gate array "A" synthesizes the amplitude modules, while the gate array "B" 
synthesizes the phase. 
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Figure 3 

Complete (upper) and simplified (lower) diagrams of states representing the purifica- 
tion of the density matrix of a single-qubit system. Starting from the input state |00), 
the active information is processed along thick lines, while thin lines correspond to 
absence of information. 



For the sake of simplicity, we consider an evolving quantum system con- 
stituted by a single qubit of density matrix p, and by a single ancillary qubit 
on which the partial trace operation is performed. This choice causes no loss 
of generality, as the Kraus representation can easily be generalized for main 
and auxiliary systems of higher dimensionality by appropriate immersions 
and permutations of diagrams of states [[1]. 

The processing of information involved in the time-evolution is illustrated 
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Figure 4 

Kraus representation for tiie time-evolution of a general density matrix. The overall 
number of qubits is preserved during the evolution process represented by the unitary 
matrix f/: 71' + m' = n + m. 



by the diagram of states in Figure |5] According to Figure |4] the overall num- 




Figure 5 

Quantum circuit (left) and diagram of states (right) representing the Kraus time- 
evolution of a single-qubit main quantum system, supported by a single-qubit ancillary 
system in the most significant position. In the diagram, the active lines of the unitary 
matrix U determine the evolution of the overall system. A similar representation can 
be obtained for the ancillary qubit in the least significant position, by simply permut- 
ing the qubits constituting the system. 



ber of qubits is preserved during the evolution process, and in the simplest 
case considered here we have n = m = n' = m! = \. The initial density 
matrix p'^" is given by: 

/" = |0)(0|®p. (12) 

The evolution of the overall density matrix p'^" is then obtained by applying 
the unitary matrix I] appropriately divided into sub-matrices of dimension 
2x2 and, subsequently, by partial tracing on the most significant qubit. The 
final density matrix p' is: 

p' ^Ttmsb{U P^'' t/^} = 
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At 

St 



ApA"^ + CpC^ 



(13) 



with the condition of trace preservation: 



A^A + C^C ^ 1. 



(14) 



3 SINGLE-QUBIT DECOHERENCE AND ERRORS 

In any reahstic physical implementation for quantum computers, quantum 
systems which process information must always be considered open, that 
is, never perfectly isolated from the environment. This gives rise to the 
well known phenomenon of decoherence, with which we here denote any 
quantum-noise process due to the unavoidable coupling of the main system to 
the environment (see, e.g., 161, pages 335-351). In addition to decoherence, 
many sources of imprecision or perturbation have to be taken into account 
in all quantum information processing tasks, since any experimental imple- 
mentation of quantum states and operations is unavoidably imperfect to some 
extent. 

To our best knowledge, error models in the literature do not offer a com- 
plete and realistic description of the most general noise-transformations for 
single-qubit systems. In ||5l we considered all physically possible single- 
qubit errors: A general transformation of a single-qubit density matrix is de- 
termined by a set of parameters, with which we associated basic transforma- 
tions, or quantum noise channels. Here we provide a detailed analysis and 
complete description of this single-qubit error model, representing all noise 
channels by means of quantum circuits, diagrams of states and geometrical 
visualizations of their action on pure quantum states. 

3.1 General transformation of a single-qubit density matrix 

It is well known that any quantum state: 



\^)^a\0)+P\l) 



a 



\a^\ + \P'\^l, 



(15) 



can be rewritten in spherical coordinates: 





(16) 
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and represented by the Cartesian coordinates of the three-dimensioniil space 
embedding the Bloch sphere. 

The quantum state 1^*) can also be associated to its density matrix p = 
I I which, in turn, can be expressed by means of Pauli operators cr: 



1 + Z X-iY 
X + iY l-Z 



X 
Y 
Z 



(17) 



where A is the vector of the Bloch sphere coordinates of the state l^"). 

According to the Kraus representation, the most general evolution of a 
single-qubit density matrix is given by: 

p' = Y.^ipPl Y.PlFi=h (18) 



which leads to an affine transformation in the Bloch sphere coordinates of the 
state I*): 

A' = MA + c, (19) 

where M is a real matrix of dimension 3x3 and c is a translation vector of 
dimension 3. Thus, twelve parameters are necessary and sufficient to charac- 
terize a generic quantum-noise operation acting on a single-qubit system. 

These twelve parameters can not be chosen arbitrarily, since they are sub- 
jected to the mathematical conditions of complete positivity of the quantum- 
noise transformations, that is, we always need to assure the physical feasi- 
biUty of the transformations determined by the parameters. This restriction 
may require a very difficult treatment if we choose to rely only on analytical 
descriptions. On the contrary, by studying the parameters by means of equiv- 
alent quantum circuits, the mathematical conditions for completely positive 
transformations are inherently and straightforwardly satisfied, and the same 
method easily allows for extensions to higher-dimensional cases. 

In order to develop descriptions by quantum circuits, we thus decompose 
the matrix M as one diagonal and two orthogonal matrices, M = 0\D Oj , 
and we associate the basic parameters with transformations of a general quan- 
tum state residing in the Bloch sphere. Three parameters correspond to ro- 
tations of the sphere around the axes x, y, z; three parameters correspond to 
displacements of the sphere along the axes x,y,z; three parameters corre- 
spond to deformations of the sphere into an ellipsoid with x, y or z as sym- 
metry axes. The last three parameters correspond to further rotations of the 
already deformed sphere, in order to obtain an ellipsoid with symmetry axis 
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along arbitrary directions. Thus, we only need to illustrate in detail the first 
nine transformations, since deformations of the Bloch sphere into an ellipsoid 
with symmetry axis along an arbitrary direction can be obtained by appropri- 
ate compositions of the first nine transformations. 

3.2 Rotations of the Bloch sphere around the coordinate axes 

Rotation errors are the simplest case of single-qubit errors. They are unitary 
transformations and can be obtained immediately as special cases of single- 
qubit unitary gates, whose diagrams of states are described in detail in [|T|. 
For the reader's convenience, the corresponding transformations of the Bloch 
sphere coordinates are: 

X' = X 



Rrr — 



Ry 



COS 2 

—i sin T 



cos ; 



cos : 



y = cos 6* r 

Z' = smeY - 



-sineZ 
cosOZ 



' sm 



s 1 



X' ^cosOX -sm9Z 
Y' = Y 

Z' = sm0X + coaeZ 



Rz 



I sm ■ 











cos : 



I sm ■ 



X' = coseX -sin9Y 
Y' = smeX + cos e Y 
Z' = Z 

(20) 

where the error parameter 9 denotes the rotation angle, that is, the intensity 
of the rotation error applied to the original single-qubit quantum state. 

3.3 Deformations of the Bloch sphere along the coordinate axes 
Deformation errors are non-unitary errors: Consequently, they require a phys- 
ically feasible representation involving a unitary operation in an extended 
quantum system. 

The deformation of the Bloch sphere into an ellipsoid centered at the ori- 
gin of the sphere with axis directed along x,y or z can be respectively im- 
plemented by means of the bit flip, bit-phase flip and phase flip channels. 
The corresponding quantum circuits are illustrated in Figure |6] (left), where 
the environment is represented by a single ancillary qubit. The processing of 
information caused by each deformation channel is clearly illustrated in the 
corresponding diagrams of states in Figure |6] (right). 

The ICraus operators of each channel can be easily derived from the corre- 
sponding diagram of states. For the bit flip channel we have: 



cos - 
2 



Fo = 



sm - 
2 



(21) 
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Figure 6 

Quantum circuits (left) and ttie corresponding diagrams of states (right) representing 
deformations of the Bloch sphere along the x, y and z axes. From top to bottom, the 
bit flip, bit-phase flip and phase flip channels. The ancillary qubit representing the 
environment is set to the initial state | * ) = cos ||0)-|-sin||l), with < ^ < tt . 



For the bit-phase flip channel we have: 



cos - 
2 



Finally, for the phase flip channel we have: 



cos - 
2 



sm - 
2 



sm - 
2 



,y. 



(22) 



(23) 



The Kraus operators of the three channels induce the following transfor- 
mations in the Bloch sphere coordinates of the initial state: 



X' = X 

Y' = cos9Y 

z' = cosez 



X' 


= cos 


9X 


Y' 


= Y 




Z' 


= cos 6 


)Z 



x' = cosex 

Y' = COS 6* F (24) 
Z' = Z 
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These transformations can be finally visualized in the Bloch sphere repre- 
sentation (Figure |7]i. 




Figure 7 

Graphic visualization of deformations of the Bloch sphere along the x, y and z axes 
by means of the bit flip channel, the bit-phase flip channel and the phase flip channel, 
respectively. The Bloch sphere (left) represents initial unperturbed states, while the 
ellipsoids (right) show how the Bloch sphere is affected by the noise channels, for 
increasing values of the error parameter (6 — 0, ^, j). 



3.4 Displacements of the Bloch sphere along the coordinate axes 

As happens for deformations, also displacement errors are non-unitary errors 
and, consequently, require a physically feasible representation involving a 
unitary operation in an extended quantum system. The most important feature 
of all displacement channels is that they must perform at the same time a 
correlated deformation, so that the resulting matrix p' is still a density matrix 
representing a physically feasible state. 

The transformations induced by displacements of the center of the Bloch 
sphere along the axes x,y or z can be implemented by means of the ampli- 
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tude damping channel and its appropriate modifications to reverse or rotate 
the displacement directiorQ- The quantum circuits for the standard and mod- 
ified amplitude damping channels are illustrated in Figures [8] and [TOfleft). 
where the environment is again represented by a single ancillary qubit. The 
processing of information caused by the different instances of the amplitude 
damping channel is clearly illustrated by the diagrams of states in Figures |9] 
and[IO](right). 



|0> 



-9/2 



9/2 



D 



|0> 



Figure 8 

Quantum circuits representing displacements of the Bloch sphere along the axis z. 
From left to right, the amplitude damping channel and its reversion to achieve the 
displacement in the positive and negative directions of the axis z, respectively. 



The Kraus operators for each instance of the amplitude damping channel 
can be easily derived from the corresponding diagram of states. For the stan- 
dard amplitude damping channel we have: 



1 

cos! 



sin^ 

o' 



For the reversed amplitude damping channel, we have: 









sin ■ 



(25) 



(26) 



To rotate the amplitude damping channel for displacements along the axes x 
and y, we apply respectively the unitary transformations: 



1 



1 ±1 

Tl 1 



1 



1 ±i 

±i 1 



(27) 



* We remark that the ampHtude damping channel actually desciibes the most general defor- 
mation errors: It performs the best possible deformation of the Bloch sphere in relation to the 
desired displacement of its center along the considered coordinate axis, as given by the high 
degree of tangency of the resulting ellipsoid to the initial Bloch sphere. 
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Figure 9 

Diagrams of states representing displacements of tiie Bloch sphere along the axis z. 
From top to bottom, from left to right: complete and simplified diagrams of states of 
the amplitude damping channel, and diagram of states of its reversion. 




Figure 10 

Quantum circuits (left) and diagrams of states (right) representing displacements of the 
Bloch sphere along the axes x and y. Circuits and diagrams are derived by rotating the 
displacement direction along which the amplitude damping channel and its reversion 
act. Rotations are obtained by applying the appropriate unitary matrices U and 17^, 
while "D" denotes the gate array corresponding to the amplitude damping channel, as 
shown in Figures[8]and|9] 



Thus we obtain the following Kraus operators for x displacements: 



Fn = - 
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1 - 


f cos 1 
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. ±(1 


— cos |) 




1 






~ 2 


-sin§ 



±(l-cos|) 



cos ; 



± sin \ 



(28) 
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and for y displacements: 



1 + COS 

T^(l 



2 

COS ' 



±z(l-cosf) 



±i sin 
sm 



2 """^ 2 
=F« sin ! 



(29) 

2 -r"^"-2 - 

The Kraus operators induce the following transformations in the Bloch 
sphere coordinates of the initial state, respectively for the negative and posi- 
tive directions of the axes x,y, z: 



X' = =F sin^ ± cos^ OX 

Y' = cos er 
z' = cos ez 



X' = cos ex 

Y' = T sin^ ± cos2 OY 

z' ^ cos ez 



X' ^ cos ex 

Y' = cos eY 

Z' sin^ e ± cos2 eZ 



(30) 



These transformations can be finally visualized in the Bloch sphere repre- 
sentation in Figure [in 

3.5 The depolarizing channel 

After describing the complete set of basic single-qubit noise operations, we 
conclude this section by mentioning a well known and widely used single- 
qubit decoherence model: the depolarizing channel (see, e.g., ||6|, pages 340- 
343). The depolarizing channel can be described by means of the quantum 
circuit in Figure [12] (left). 

The processing of information performed by the depolarizing channel is 
clearly illustrated in the diagram of states in Figure [12] (right). The main 
single-qubit system is coupled with the environment, which is now repre- 
sented by two ancillary qubits, set in the pure state: 



1*) = a|00) +/3|01) +7|10> + 



with the normalization condition: 



\a\' + m' + \l\' + \S\' = l. 



(31) 



(32) 



The channel applies one of the the operators I,ax,(yyOr to the least signif- 
icant qubit if the two most significant qubits are in one (or in a superposition) 
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Figure 11 

Graphic visualization of displacements of the Bloch sphere along the x, y and z 
axis. The Bloch sphere (left) represents initial unperturbed states, while the ellipsoids 
(right) show how the Bloch sphere is affected by the noise channel, for increasing 
values of the error parameter (9 = 0, ^, j, 
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Figure 12 

Quantum circuit (left) and diagram of states (right) representing a generalized depo- 
larizing channel. Depending on the choice of the parameters \a\, |/3| and |7|, several 
noise operations can be described as special cases of this channel. 



of the states |00), |01), |10) or |11). The final state of the single-qubit system 
is obtained by tracing over environmental qubits. 

Depending on the choice of the parameters |a|, \f3\ and I7I, several noise 
operations can be described as special cases of this channel (for instance, the 
bit flip or the phase flip channels). 

The standard depolarizing channel is obtained when: 

\a\^ = cosH, = = |5|2 = (33) 
In this case, the Kraus operators become: 



^/3 



CTt, 



i.(3) = !j^.,, F(^) = ^.. (34) 

and the Bloch sphere coordinates of the initial state are transformed as fol- 
lows: 

X' = (1 - I sin2 0) X 

F' = (1 - I sin^ e) Y (35) 
Z' = sin^ e) Z 

Thus, the Bloch sphere is deformed into an elUpsoid, centered at the origin of 
the Bloch sphere, whose axes are directed along x, y and z, the deformation 
rate being the same along each axis. 



18 



Thus, the depolarizing channel can not be considered as a general model to 
describe and study single-qubit decoherence effects, since it only includes a 
small subset of specific noise operations, although involving an higher dimen- 
sional auxiliary system in respect to the previously illustrated single-qubit 
noise channels. 

4 CONCLUSIONS AND FUTURE DEVELOPMENTS 

We have explored the main processes involved in evolution of general quan- 
tum systems by means of Diagrams of States, a novel method to graphically 
represent and analyze how quantum information is elaborated during com- 
putations performed by quantum circuits. We have offered complete and de- 
tailed descriptions by diagrams of states of partial trace operations, density 
matrix purification and evolution by Kraus operators, single-qubit decoher- 
ence and noise channels. 

In our opinion, diagrams of states can be used as an auxiliary or as an 
alternative approach to standard methods, both to investigate and to conceive 
quantum computations. Analytical study and Feynman diagrams alone are 
often too synthetic to clearly visualize how quantum information is processed 
by computations. On the contrary, the dimension of the graphic representation 
of states grows exponentially in respect to the dimension of the examined 
quantum system, thus offering a complete and detailed visualization of the 
computational process. 

Diagrams of states appear to be most useful whenever the quantum opera- 
tions to be analyzed are described by very sparse matrices, since only non-null 
entries of matrices are associated with diagram lines which contain active in- 
formation. This way, the resulting diagrams show clearly and immediately 
the significant pattern along which quantum information is processed by the 
computation from input to output. Indeed, several quantum computations ac- 
tually involve operations satisfying this requirement, and evidence of this is 
also provided by the processes illustrated in this paper 

Further computations are going to be explored by this graphic representa- 
tion, among which quantum algorithms |,9J and entanglement-based practical 
applications ifTOl . 
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